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ABSTRACT 

We derived one-fluid equations based on a relativistic two-fluid approxima- 
tion of e 1 * 1 pair plasma and electron-ion plasma to reveal the specific relativistic 
nature of their behavior. Assuming simple condition on the relativistic one-fluid 
equations, we propose generalized relativistic magnetohydrodynamic (RMHD) 
equations which satisfy causality. We show the linear analyses of these equations 
regarding various plasma waves to show the validity of the generalized RMHD 
equations derived here and to reveal the distinct properties of the pair plasma 
and electron-ion plasma. The distinct properties relate to (i) the inertia effect 
of electric charge, (ii) the momentum of electric current, (iii) the relativistic 
Hall effect, (iv) the thermal electromotive force, and (v) the thermalized energy 
exchange between the two fluids. Using the generalized RMHD equations, we 
also clarify the condition that we can use standard RMHD equations and that 
we need the distinct RMHD equations of pair and electron-ion plasmas. The 
standard RMHD is available only when the relative velocity of the two fluids is 
nonrelativistic, a difference of the enthalpy densities of the two fluids is much 
smaller than the total enthalpy density, and the above distinct properties of the 
pair/electron-ion plasma are negligible. We discuss a general relativistic version 
of the equations applicable to the pair and electron-ion plasmas in black hole 
magnetospheres. We find the effective resistivity due to shear of frame dragging 
around a rotating black hole. 

Subject headings: plasmas, relativity, methods: analytical, galaxies: active, galax- 
ies: jets, galaxies: magnetic fields, galaxies: nuclei 
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Introduction 



The most powerful engines in Universe are located in centers of active galactic nuclei 
(AGNs), quasars (QSOs), micro quasars (/xQSOs ), and gamma-ray bursts (GRBs). Observa- 
tions have shown all of them eject relativistic je t s (IPearson fc Zensus 1 987; Bi retta et al. 1999 
Junor et al. 19991 : iMirabel fc Rodriguez 1994 ; iMirabel fc Rodriguez 19981 : iKulkarni 1999h . 



In spite of the drastic difference of their characteristic scales and powers, it is believed that 
the activities of these objects, such as the relativistic jet ejecti on, are supported commonl y 
by the drastic phenomena of accretion disks around black holes ( IMirabel fc Rodriguez 19981 ). 
However, the distinct mechanism of the activity is not confirmed yet. Models with inter- 
action between the plasma and magnetic field in the very strong gravity of the black hole 
are thought to be most promising. The region where the plasma and magnetic field interact 
each other near a black hole is called a black hole magnetosphere. A black hole magneto- 
sphere consists of an accretion disk and a corona around a black hole and it is thought to 
be composed of various kinds of plasmas. In a cas e of an AGN, it is suggested that the disk 
is made of electron-ion plasma (IFord et al. 19941). the core of the relativistic jet is mainly 
of electron-positron (pair) plasma (j War die et al. 19981 ). and the corona is by both of them, 
while the actual components of such plasmas have not been confirmed observationally yet. 
It is a natural idea that a corona near a black hole consists mainly of pair plasma, and 
near its disk, mainly of electron-ion plasmas. Such a difference of the kinds of plasmas may 
influence the dynamics of a black hole magnetosphere, but none of the works has been done 
to investigate it. 

Black hole magnetosphere has been studied on the basis of the non- relativistic m agne- 
tohydrodynamics (MHD) with pseudo-Newtonia n potential (IPaczvnski &: Wiit a 1980) an d 



Takahashi et al. 1990 



Tomimatsu 19941 ) . 



ideal general relativistic MHD (GRMHD) (e.g. 
The relativistic MHD (RMHD) or GRMHD is a one-fluid approximation of the plasma 
and is based on the relativistic conservation laws of particle number, momentum, and 
energy, Maxwell equations, and the simple Ohm's law in the plasma rest frame. The 
simple Ohm's law means that the electric field E' is proportional to the current den- 
sity J' in the plasma rest frame: E' = rjj' where rj is the electric resistivity. We call 
such RMHD (GRMHD) equations the "standard" RMHD (GRMHD) equations in this pa- 
per. When we set rj zero, the standard RMHD (GRMHD) becomes the "ideal" RMHD 
(GRMHD). Recently, numerical simulations of the ideal GRMHD have become prevailing 
and yielded i nteresting, impor tant results with respect to relativistic outflow/jet formation 
(IKoide 20041: iMcKinney 20061) and e xtraction of black hole rotational energy by magnetic 
field ( IKoide et al. 20021 : IKoide 20031 : iKomissarov 20051 ). However, since they used the stan- 
dard GRMHD and thus the one-fluid approximation, they ignored the unique properties of 
the pair/electron-ion plasma. To begin with, we have to check the validity of the standard 
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RMHD and clarify its application conditions for plasmas around black holes. 

To clarify these points, we reconsider the RMHD equations. Such a task was first 
performed by Ardavan (1976) using the Vlasov-Boltzmann equation for a pulsar magneto- 
sphere. It yielded a relativistic version of the generalized Ohm's law and a new condition 
for the validity of the MHD approximation for a pulsar magnetosphere (where Lorentz fac- 
tor of plasma is much larger than unity). A more generalized treatment, which included 
annihilation of electrons and positrons, radiation, Compton scattering, and pair photopro- 
duction was formulated by Blackman & Field (1993) and Gedalin (1996). Reconsideration 
of ideal MHD in a neutral cold plasma based on two-fluid approximation was presented by 
Melatos & Melrose (1996), who investigated the conditions under which the MHD approxi- 
mation breaks down. To investigate black hole magnetospheres, Khanna (1998) formulated 
a general relativistic version of the two-fluid approximation in the Kerr metric. A more gen- 
eralized version in a time-varying space-time was derived by Meier (2004) from the general 
relativistic Vlasov-Boltzmann equation. 

To investigate the distinct properties of the pair and electron-ion plasmas and the appli- 
cability of the standard RMHD, we derived one-fluid equations using a relativistic two-fluid 
approximation of a plasma consisting of positively charged particles, each has charge e and 
mass m + , and the other fluid of negative particles, each with charge — e and mass m_. 
Throughout this paper, we assume the charge of ions is e, while generalization for plasmas 
with arbitrary charge of ion is not difficult. In an electron-ion plasma case, m + is much 
larger than m_, while in a pair plasma case, they are equal. With respect to the one-fluid 
equations with finite resistiv ity larger than a certain value, causality seems to be broken 



mathematically (IKoide 20081 ). We consider the dispersion relation of the electromagnetic 



waves and show the mathematical condition where the equations satisfy causality. 

On the basis of the above results, we propose the "generalized RMHD equations" for 
the plasmas, which are introduced in this field for the first time by the present work. These 
equations can reveal the distinct properties of the electron-ion and pair plasmas. The differ- 
ences between the generalized RMHD equations and the standard RMHD equations concern 
(i) the inertia effect of charge, (ii) the momentum of current, (iii) the relativistic Hall effect, 
(iv) the thermal electromotive force, and (v) the thermalized energy exchange rate. Mainly 
using linear analyses, we investigate the influence of such distinct properties to plasma dy- 
namics which cannot be clarified from the standard RMHD equations. The conditions that 
the distinct properties of the electron-ion and pair plasmas are negligible and the premise 
conditions of the generalized RMHD clarify the application condition of the standard RMHD 
equations. 



In section [2J we derive one-fluid equations from the relativistic two-fluid equations. In 
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section [3], dispersion relation of electromagnetic wave is shown on the basis of the one-fluid 
equations and we confirm the mathematical condition of causality of the one-fluid equations. 
In section HI we propose a set of generalized RM HD equations with any mass ratio m + / m_ 
from the one-fluid equations. To show the reasonable property of the generalized RMHD 
equations, in section [5j we investigate the linear modes of various plasma waves derived by 
the generalized RMHD equations and show the unique nature of the pair and electron-ion 
plasmas. Section [6] shows the unique properties of the electron-ion and pair plasmas and 
estimations of them. In section [7J we discuss the break-down condition of the generalized 
RMHD equations and the non-linear distinct nature of the electron-ion and pair plasmas in 
the jet forming region around the black hole. In section [HJ we suggest the distinct general 
relativistic effects shown by the generalized GRMHD equations. Here, we indicate that shear 
of frame dragging induces effective resistivity near a rotating black hole. In section the 
results are summarized. 



We employ relativistic two-fluid equations of plasma consisting of two fluids. One fluid 
is composed of positively charged particles, and the other is composed of the negatively 
charged particles. We derive a set of one-fluid equations equivalent to the two-fluid equations 
in the Minkowski space-time (x°, x 1 , x 2 , x 3 ) = (t,x,y,z) where the line element is given by 
ds 2 = —(dx ) 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 = r] fJiU dx fJ, dx 1 ' . Throughout this paper, we use units 
in which the speed of light, the dielectric constant, and the magnetic permeability in vacuum 
all are unity: c = 1, e = l, fio = 1. The relativistic equations of the two fluids and the 
Maxwell equations are 



where variables with subscripts, plus (+) and minus (-), are those of the fluid of positively 
charged particles and of the fluid of negative particles, respectively, n± is the proper particle 
number density, U± is the four-velocity, p± is the proper pressure, h± is the relativistic 
enthalpy density, is the electromagnetic field tensor, *F^ U is the dual tensor density of 
F^, i? M is the frictional four-force density between the two fluids, and J M is the four-current 
density. We will often write a set of the spacial components of the four-vector using a bold 



italic font, e.g., U± = J = (J 1 , J 2 , J 3 ), R = (R\R 2 } R 3 ). We further define 



2. Two-fluid equations and one-fluid equations 



d v {n ± U v ± ) 
d v {h ± U$U v ± ) 
8 U *F" V 
d v F» v 



0. 

J", 



o, 

-«9>± ± en±ri^U^F au ± iT, 



(1) 
(2) 
(3) 
(4) 
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the Lorentz factor 7± = U±, the three- velocity V± = U±/^±, the electric field = F°\ the 
magnetic flux density Bi = Ylijk e ijkF-' k /2 (e^ is the Levi-Civita tensor), and the electric 
charge density p c = J°. Here, the alphabetic index (i,j,k) runs from 1 to 3. Throughout 
this paper, we assume that the two fluids are heated only by Ohmic heating and disregard 
nuclear reactions and pair creation and annihilation. We also ignore radiation and quantum 
effects. 

To derive one-fluid equations of the plasma, we define the average and difference vari- 
ables as follows: 



P 

n 



m + n + + m_n_, 
P_ 

7X1 
1 



17" = -(7xi + n + U» + m-n-U»), 



(5) 
(6) 

(7) 
(8) 



jn = ein+Ut-n-U^), 
where m = m + + ra_. Then we can write 

n±m = -( P U^± ^ jA . (9) 

Txi V e / 

We also define the average and difference variables with respect to the enthalpy density as 

h± h 



h) 
A/i 1 " 

h* 
Ah* 
Ah} 



n 



= n 



nt nt 



n 



= n 



= n 



■h + 


2m_ 


h- 2m + 


n\ 


m 


n 2 m 


h+ 




1 nt 


n\ 


V m / 


'K 


2m + 


h- 2m_ 


n\ 


m 


n 2 m 


K 


/2m_ N 


)-- 


nt 


V m , 


1 nt 



m 



2m_ 



m 



(10) 
(11) 
(12) 

(13) 
(14) 



We find the following relations between the variables with respect to the enthalpy density, 



h x = 4/i/i 1 " - A/j,Ah\ 
Ah* = 2A/i/i t + A/i 1 ", 
Ah} = SfiAfih* + 4(1 - 3fi) Ah\ 



(15) 
(16) 
(17) 
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where \x = m + m^/m 2 and Afj, = (m+ — m_)/m. It is noted that we have the relation 
H — [1 — (A/i) 2 ]/4. We also find the relations with respect to the average and difference 
variables of the two fluids as 



m + n + U t t + m_n_C/^ = pU^ (m + — m) J M , 

e 



h+UtUl + h-UtVL 



hX A At 

h ] u»u v + — — —j^j v + ^—{u^r 



(2ne) 



2ne 



(18) 
(19) 



hi 



■p(U ti J v + J»U h 



A/ ^ P 2 fw^ + ^A/^(jm) 



4n 2 e Iran 2 
Using the above variables and their relations, we derive the one-fluid equations from the 
two- fluid equations (pQ), <^ as 



d, 



0. 



/,t A/jt 



(2ne) J 



2ne 



(21) 

-d^p + J u F^ v , (22) 



2ne 



At A/jfl 

' (17" J* + .W) + AtfUW + ^—^J^F 



2ne 
1 



(2ne) 



2ne 

where p = p + + p- , Ap = p+ — p- , and 



Ap 



0"(A/xp -Ap)+[U U -J v F^ v + — 



n e 



M y 1 



(23) 



(24) 



Here, M = U U U U , Q = U U J U , and G is the thermal energy exchange rate from negatively 
charged fluid to the positive fluid (see Appendix [SJ . In the electron- ion plasma case, we 
usually assume that the energy exchange due to the friction between the two fluids is not 
significant: G = 0. On the other hand, in the pair plasma case, since the energy of electron 
fluid and positron fluid is exchanged frequently, we set 

2m 2 Q 2 - MW 



G 



(2mnM) 2 - (mQ/e) 2 



0. 



(25) 



where W = J U J U and 9 is the redistribution coefficient of the thermalized energy due to the 
friction of the two fluids (see Appendix lAl. 



3. Propagation of electromagnetic waves and causality 



With respect to the one-fluid equations of a pair plasma, Koide (2008) indicated that 
causality with the electromagnetic wave propagation breaks down if the electric resistivity 
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of the pair plasma exceeds a critical value. To clarify the causal condition on the one-fluid 
equations, we investigate electromagnetic waves in unmagnetized, uniform plasma using the 
linear analysis of the one-fluid equations f[2"Tj) - ([2"3"j) and the Maxwell equations ([3]) and (TjJ. 
We write the variables of the uniform, rest plasma and uniform magnetic field as p = p, 
p = p, h> = h, U = 0, B = 0, and E = 0. When the perturbations of the variables are 
written by p = p — p, p = p — p, h = h* — h, U = U, B = B, and E = E, we have the 
linearized equations in the three-vector form, 

jLp + pV-L/ = 0, (26) 

§- t (h~P) = -hV-U, (27) 

h^U = -Vp, (28) 

h d X 

p--K = — V ■ (App - Ap) + E - 7) J, (29) 

q at 2q 

V-E = p e , (30) 

V B = 0, (31) 

|b = -VxE, (32) 

J+^-E = VxB, (33) 
at 

where q = ne is the "latent" charge density, K = J/(ne), q = ne, and K = J /(fie). 
Hereafter, we note an equilibrium variable by a bar, and perturbation by a tilde. We assume 
the resistivity is uniform and constant. We also assume perturbation of any variables A is 
proportional to exp(ik-r — iut) = exp(ip AtI/ /c /i x !/ ), where = (u, k) is the constant covariant 
vector called the wave number four-vector. Furthermore, for simplicity, we investigate the 
transverse modes of the electromagnetic waves in non-magnetized plasma, thus we set 

E, B, U ± k. 

Then we have following linearized equations, 

p = , ~h-p = 0, (34) 

-iuMJ = 0, (35) 
h ~ 

—iuop—J = E — rjJ, (36) 
q 2 

-iuB = -ik x E, (37) 

J - iluE = ik x B. (38) 
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These equations yield 



77 - iw/i^J J = E, (39) 
iujJ + u 2 E = k 2 E. (40) 



Then we obtain the dispersion relation 



7) — iu—u;] (k 2 — lu 2 ) = iuj. (41) 
T J 

When we write k = k/{r]u 2 p ) and u = uj/(r]u 2 p ) using u> 2 = q 2 /(fih), we have 

H{k 2 -Cj 2 )(1-iuj) =iu, (42) 

where H = {j]uj xp ) 2 = (r/q) 2 / (fih). This equation is mathematical ly identical t o the dispersion 
relation of electromagnetic wave in the resistive pair plasma (IKoide 20081 ). The detailed 
analysis of the dispersion relation (1421) performed by Koide (2008) shows that the group 
velocity of the electromagnetic wave is smaller than the speed of light when H < 2, while it 
is larger than the light speed (unphysical) when H > 3. That is, if > 3 is forbidden, while 
the details in the range 2 < H < 3 are shown in Koide (2008). Consequently, the condition 
of the causal one-fluid equations is given by H < 2, that is, 



77 < — . (43) 



LO rp 



4. Generalized RMHD equations for electron-ion and pair plasmas 



For simplicity, we set the assumptions, 

A/i f = 0, M — U U U U = -1. 
When we write h = h) , we find 

h x = 4/j/i, Ah x = 2Afih, A/i" = SfiAfih. 
From equations (f2"Ti) - (l2"3"|) and the above assumptions, we obtain 
d v {pU") = 0, 



h ( u^u u + 4 w 



Q 



^(U^F + .FU h 



= -<9> + J U F» V1 
2fiAfih 



-.FJ l 



q 



^(Afip - Ap) + [U v - ) F\ - V [J^ + Q(l + e)u" 

2q \ Q 



(44) 

(45) 

(46) 
(47) 



(48) 
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These equations with Maxwell equations ([3]) and (j4j) can be regarded to satisfy causality when 
H = (rjq) 2 1 '(fih) < 2 everywhere. We call equations fj46|) - fj48|) with H < 1 the "generalized 
RMHD equations." 

The generalized RMHD equations for the electron-ion plasma are given by the limit 
m « m + ^> m_ and an assumption 9 = as 



mq 



h 



d v {pU v ) 
d u {hU»U v ) 



o. 



-Vp- + (V - -J u ) F» v - rj[J- 



These equations (|49p - (|5ip are similar to the standard RMHD equations 

d v {pU v ) = 0, 
d u {hU»U v ) = -d fM p+ J U F» V1 
U U F^ U = ri[J^ + QU"]. 



(49) 
(50) 

QU"]. (51) 



(52) 
(53) 
(54) 



Especially, the equations with respect to the mass density p (equations (|4"9"1) and (1521 ) and 
the momentum density (equations (1501) and (|53|) ) are identical except for the second term 
of the left-hand side of equation (T47j). The difference between the two sets of equations are 
found only in the Ohm's laws (j5ip and (l54"j) . In the standard RMHD equations, inertia of 
current density, thermal electromotive force, and the Hall effect are ignored. 



The generalized RMHD equations for a pair plasma are given by setting m + = m_ 
(m e is the mass of an electron) as 



d v { P U v 



0. 



h (uw + 



(2?) 



-d^p + rF» v , 



= m e 

(55) 
(56) 



Q 



— d^Ap + U u F^ u -r][J" + Q{l + Q)U^}. (57) 
2q 



In this pair plasma case, the Hall term disappears, while the current inertia effect is found 
in the equation of motion (1561 . The thermal electromotive force is different from that of 
electron-ion plasma. It is also noted that G is not negligible. 



In this section, we proposed a set of generalized RMHD equations (l4"Bl - (j4"81) . the electron- 
ion plasma RMHD equations (14"9"]) - (15T|) . and the pair plasma RMHD equations (I55|) - (157|) . 
The comparison between these sets of equations and the standard RMHD equations (152]) - 
( I54p reveals the unique properties of the electron-ion and pair plasmas. We examine these 
unique properties in the following sections. 
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5. Linear analyses of RMHD waves 

We newly proposed a set of generalized RMHD equations (}4"6"|) - (j4"8"}) which are applicable 
to electron-ion and pair plasmas in the previous section. In this section, we show the linear 
analyses of these equations concerning various plasma waves to show the reasonable property 
of the generalized RMHD equations and to clarify the distinct properties of the electron- 
ion and pair plasmas. First, we derive the general dispersion relations with respect to the 
generalized RMHD equations, and next we reveal the unique nature of the electron-ion and 
pair plasmas, respectively. 

We investigate waves propagating in a uniform, rest plasma and a uniform magnetic 
field. Linearized equations of perturbations, p = p — p, p = p — p, h = h/ — h, U = U, 
B = B — B, and E = E, are given by 

^p + pV-C/ = 0, (58) 

- d ~ 

h—U = -Vp + JxB, (59) 
at 

p——K = — V ■ (App - Ap) + (U - ApK) xB + E- V J, (60) 
q at 2q 

V-E = p e , (61) 

VB = 0, (62) 

= — V x E, (63) 

J + ?-E = VxB. (64) 
at 

These equations are closed with the equation of state (EoS), h = h(p,p) = ph(p/p). The 
adiabatic EoS for single- component relativistic fluids, which are in thermal equilibrium, has 
been known, and is given by 

h± K 3 (p±/p±) _~ 
p± K 2 {p± p±) 



(IChandrasekhar 19381 ; ISynge 19571 ). Here, K 2 and K 3 are the modified Bessel functions of 



the second kind of order two and three, respectively. Useful approximation of the EoS 

3p + 3 
3p + 2 



h RC (p) = l+p(l + ^±^-) (60) 



was proposed by Ryu et al. (2006), where p = p/p. This expression has the functional 
dependence similar to h(p) of equation fl65|) and replaces complex calculations of modified 
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Bessel functions by algebraic calculations. When we consider the adiabatic one-component 
fluid in the rest frame, equations ([I]) and (T5]) reduce to 

d v {n ± U v ± ) = 0, (67) 
d v {h±U$U v ± ) = -d»p±. (68) 

The linearized equations of the perturbations, n± = n± — n±, p± = p± — p±, h± = h± — h±, 
and U£ = U£- U£ where U£ = (1, 0, 0, 0) are 

d n± + n±^d i U i ± =0, (69) 

i 

d ~h ± + h ± Y,diU l ± =-d° P± , (70) 
h±d Ui = -d* P± . (71) 

These equations yield 

h+ h± — p± , . 

— = - . 72 
n± h± 

From equation ( |65i) . we find 

^= h '^ =r(P±). (73) 

P±P± i-h'(£±\ \p±J 

In general, V is not constant and is a function of p±/p±. V is called the effective adiabatic 
indexjj] When T (p + /p + ) ta T(p_/p_) which is satisfied only if p± <C p± or p± p± (e.g., 



Rvu et al. 20061 ) J 2 I we have (1/T)p = Tp, that is, 



V P V \dp J \p 
In general, (p/p)(p/p) = T(p,p) depends on both p and p. 



PP Pfdp\ T (P\ (74) 



5.1. Plasma oscillation 

To show the inclusion of plasma oscillation as a fundamental mode of the generalized 
RMHD equations (H61) - (ll8l with the Maxwell equations (jBJ), (BJ, we simply assume 

^ = 0,5 = 0,^ = 0. (75) 



1 The polytropic index is given by N — (T — 1) . 

2 We can check this using the approximation of /irc in equation (|66p easily. 
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Because p — 0, according to equation fl73|) . we have p = Ap = 0. The linearized equations 
fl5ED-([64D become 

V—fr = E\ (76) 
q at 

X + ^-E 1 = 0. (77) 
at 

Immediately, we have the dispersion relation 

u 2 = u 2 p , (78) 

where u 2 p = q 2 / (pih). This mode is recognized as the plasma oscillation, because u rp is iden- 
tical to the plasma frequency u p = \Jne 2 / {Jj/m) when p = 0: for the electron-ion plasma, the 
plasma frequency is oo ip = ^/ne 2 /m e , and in the pair plasma case, ua ip = y^2ne 2 /m c . As we 
will show in next subsection, u rp is the plasma frequency when the pressure is relativistically 
high. 



5.2. Compressional modes 

In this subsection, we derive a dispersion relation of a longitudinal oscillation modes 
(U || k) in an unmagnetized, rest plasma with uniform, finite pressure p. For simplicity, we 
assume the temperatures of the two fluids are the same: T = T + = T_. Using equation (jSJ) 
and 0-th component of equation (jSJ), we have •yn± = (jp ± m T p c /e)/m when 7 = 7+ « 7_. 
Using these equations, we have 



d'jn 



-— V • (pC7) + — (V • J). 

m e 



(79) 



In the present non-relativistic case, setting 7 = 1, we have 

| = _l v .(pC/ ) + V-J). 

at m e 
If T(p + /p + ) = r(p_/p_) = r is uniform and constant, we obtain 



dp 


2rr 


di 


m 


9 A ~ 









V-(pC/)-^(V-J) 



(81) 
(82) 
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We have the linearized equations, 



q 2 dt 



-Vp, 
1 

2~q 
0. 



V(App- Ap) + E, 



These equations yield 

ph d 2 y 

q 2 dt 2 

- d 2 ~ 



_MI V (V • U) + g(l + (A,,) 2 )V(V J) -J, 
2r ^V ( V.t/)-^V(V.J). 



m e 
Because J \\ k and U \\ k in the longitudinal modes, we obtain 



(83) 
(84) 
(85) 

(86) 
(87) 



2 ph ~ 

~ w — J \\ 
q 2 



^^-I|(l + (A^) 2 )fc 2 J||-J||, 



-u 2 hU l{ 



2TTp 



m 



19) 



where J\\ = (J ■ k)/k and C/|| = (U ■ k)/k, k ^ 0. Then we get the following dispersion 
relation, 



u 2 - -HI + (A/i) 2 )A; 2 - lu 



rp 



(c, 2 -2c 2 A; 2 ) = ^(c s A;) 4 , 



(90) 



where c 2 = = (p/h)(dp/ dp). 



In a case of an electron-ion plasma (Ap pa 1, p — m c /m <C 1), when u ^> c s k/^/~p, the 
dispersion relation becomes 

(91) 



-s ;„2 



2/T ' " rp 

This expression shows the dispersion relation of the plasma oscillation for the plasma with 
finite pressure. In the case that uo 2 <C uj 2 / p, we have the dispersion relation of sound waves 



u 2 = 2c 2 s k 2 . 

In a case of a pair plasma (Ap = 0, p = 1/4), we have two modes 



2 2 i o 2 ; 2 

uj =u rp + 2c s k , 



and 



a/ = 2cik 



,2;„2 



(92) 

(93) 
(94) 
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The former is the dispersion relation of plasma oscillation and the latter is that of the sound 
wave. 

We investigate the compressional wave of the magnetized plasma. When k || B and 
U || B, the dispersion relation is the same as that of the non- magnetized plasma wave. 
Then, we investigate the case that k J_ B, k J_ B, k \\ U, and p c = in this paragraph. 
Because p e = 0, we have Ap = 0. In the case that k <C u rp and to <C u rc , the left-hand side, 
first term and the Hall term of the right-hand side of the Ohm's law (1601) are negligible. The 
validity of this assumption is checked bellow. Here, we defined a characteristic frequency of 
the Larmor rotation of an averaged particle of the two-fluid model by 



eB p 
m h 



(95) 



If we set the pressure to be zero, u rc reduces to the cyclotron frequency of the charged 
particle with mass m and charge e in the magnetic field B, u c = eB/m. Using linearized 
equations (|58l) - fl64l) . we have the dispersion relation, 



to 



2 ;„2 



vfk 



(96) 



where v\ = (Tp + B 2 )/(h + B 2 ). This mode shows the dispersion relation of the fast wave. 
It is also noted that Vf < 1. Below, we check the validity of the assumption k <ti uj rv and 
uj oj vc that we used. Using the Ohm's law (!60l) with this assumption and the dispersion 
relation (19"6"|) . the left-hand side of the Ohm's law (I60p becomes 



lui—J 



k 



rp 



[vj - l)U x B. 



(97) 



This clearly shows that the term of the left-hand side can be ignored in equation fl60l) when 
k <ti uj rp . The absolute value of the summation of the thermal electromotive force term and 
the Hall term of the Ohm's law (1601) is estimated by 



—ikAap 
2q 

A ^ T~T 
Q 



ApK x B 



Q 



5 + B z I ^ - B 



1 



B 



(98) 
(99) 



This is negligible compared to | U x B \ 
when to ^ uj tc . 



UB because the ratio of the two terms is ~ Apu/uj r 
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5.3. Transverse wave propagating along magnetic field 



We investigate transverse waves propagating through the ideal MHD plasma along the 
magnetic field lines, 

B || k, E, B, U || k, r] = 0. (100) 

We assume any perturbation A is proportional to exp(ik ■ r)A{t). The linearized equations 
become 



fih d ~ 
~Jdi J 
ik ■ E 
ik B 



= J x B, 

= {qll - AfiJ) x B + qE, 

= o, 

= o, 

= -ik x E, 

= ik x B. 



From these linearized equations, we have 
fih d 2 



d 2 ~ 



where □ = d v d v . When we set 



dJ_ 

~dt 



iujj + Q x J, 



(101) 

(102) 

(103) 
(104) 

(105) 
(106) 

(107) 
(108) 



we obtain the following dispersion relations, 



S^k' 2 = qu >2 + ^(cu' 2 k f2 -Acu 2 n 2 )-Afi(k' 2 -2uj 2 )n ]] B, (109) 
h q 



2fin 



21 



_-(k' 2 -^)+q + 
q h 



A/i(k' 2 - 2Vt 2 )B, 



(110) 



where k' 2 = k 2 — uj 2 — Vt 2 and to' 2 = cu 2 + f2 2 



Let us consider the case of pair plasma (fJ> — 1/4, A/i = 0). In this case, we have £1 = 0, 
uj' 2 = u, k' 2 = k 2 — u 2 , and the dispersion relation becomes 



CO 



B 2 

-T+ , 

ft V uj 



rp 



(A; 2 - cu 2 ) = ^ 



(111) 
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When u 2 p 3> oj 2 , we have 



u 



UJ 



A ;„2 



1 + ul ' 



(112) 



where u\ = \J B 2 /h is the Alfven four- velocity. This shows the dispersion relation of the 
relativistic Alfven wave. 

For the electron-ion plasma case (// ~ m e /m <C 1, A/i ~ 1), we arrange the dispersion 
equations fTT09|) . fTTTOj) as 



UJ. 



rp 



2(i + Ofi,i = (k 2 -2n 2 )^f. 



rp 



It is noted that there is the relation between uj tc and u Tp , 

^rp Q ' 

When fi < k, k -C w rp , ^rc *C w r p> then fin <C fc and A;' 2 = fc 2 , thus we have 



fin 



u; rc k 2 
2^1+^1' 



Substituting this equation to equation ( 1113ft . we obtain 



1 + - 



2 / ; 2 2\ 2 

When uj ic <C u> rp , k <C u> rp , we finally have 



3 ( ui rc k 



2 v<; i+< 



1 fUnp 

4 V< 



^ 2 = <(A; 2 -u; 2 ) 



(113) 
(114) 



(115) 



(116) 



117) 



lis) 



Equations ( 1116ft . ( 1118ft yield the dispersion relation of the Alfven wave in the electron-ion 
plasma case, 



u 



UJ 



fill 



A k 2 



1 + ul ' 

UJ TC 

2(l + <V r 2 p 



k 



(119) 
(120) 



It is noted that the polarization of the Alfven wave rotates, while the rotation angular 
velocity fiy is very small. 
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We examine the electromagnetic waves propagating along the magnetic field lines. From 
equations ()109p . (IllOp . we have 

uj' 2 k' 2 - Au 2 VL 2 + uj 2 p uo' 2 - ^uj 2 c k' 2 = A/juj^k' 2 -2w 2 )% (121) 
2n l] [k' 2 - to' 2 + cu 2 p + fxco 2 c ] = A^LU TC (k' 2 -2fi 2 ). (122) 

When u} ^> u rp ,uj rc , the medium of the electromagnetic wave becomes nearly vacuum, and 
we assume 

u > k' > n. (123) 



Then using u' = u, we have 



u; 2 (fc' 2 +u; r 2 p ) = (A/i)V c fc' 2 , (124) 



Tp7 

-2n {l uj 2 = Afiuj vc k' 2 . (125) 

Thus, we have the dispersion relations, 

u 2 = k 2 + u%, (126) 

fill = Aai-^H. (127) 

The finite fin shows the Faraday rotation of the electromagnetic wave along the magnetic 
field line in the electron-ion plasma. In the pair plasma case (A/i = 0), no rotation of 
polarization is found. 



6. Unique nature of electron-ion and pair plasmas 

The differences between the standard RMHD equations and generalized RMHD equa- 
tions of electron-ion plasmas and pair plasmas are only following five points. To find out 
these points, we should just consider the terms with /x and A/x in equations (p^ - lflB]) . 

1. Second term of left-hand side of equation of motion fj47|) : This term is due to the 
momentum of current. In the pair plasma case, it is significant, while it is negligible 
in an electron-ion plasma. 

2. The time-derivative term of the Ohm's law (1181) : This is due to the inertia effect of the 
current. We can ignore it only when the electric resistivity is zero in an electron-ion 
plasma. When the resistivity is finite, we cannot ignore it because of causality. 

3. The second term with respect to J M of the right-hand side of the Ohm's law (j!8|) : This 
is the Hall term. In the pair plasma case, it disappears. 
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4. The term with respect to pressure gradient of the right-hand side of the Ohm's law 
f|48j) : This is the term of the thermal electromotive force. In a pair plasma, it influences 
through Ap, while in an electron- ion plasma through p_. 



5. The term proportional to in the Ohm's law (1481) : This is due to the difference of the 
thermalized energy exchange between the electron-positron pair and electron-ion pair 
during the friction. 



The item 3 was indicated by many authors (e.g., IBlackman fc Field 19931 ; iMeier 2004 ) and 
the other items are explicitly described in this paper for the first time. As shown above, it is 
clear that most of differences come from the Ohm's law (j4"8l . The exception is the effect of 
current moment in the equation of motion shown in the distinct property [TJ The influence 
of the distinct property [T] is estimated by comparison between v / /JJ At /(ne) and U^, which 
are on the order of the square root of two terms in the left-hand side of equation (f471) . This 
implicitly contains two ratios, 

J° |J| 

J~u — : 1, J~u — : \U\ — vfnrn\U + — U-\ : \m + n + U + + m_n_L/_|. (128) 

ne ne 

To evaluate the latter ratio, we introduce the transformation between the two frames (t, r) 
and (£', r') whose relative velocity is Vq, like the Galilean transformation, 

£' = £, r b = r- v t. (129) 

For any three- vector A, the vector length \A b \ is equal to or smaller than the length of the 
vector transformed by the Lorentz transformation, \A'\. This is because 



\A'\ = \A\ { + A' ± | = | 7o (A|| - v t) + A ± \ = J|7o(A|| - ^)| 2 + |AJ 2 



> ^\A\\ - v t\ 2 + \A ± \ 2 = |A|, - v t + A ± \ = \A\ + A\\ = \A\ 

where the vector with the subscript || (_L) is its parallel (perpendicular) component to v 
and 70 = (1 — Vq)^ 1 ^ 2 . When we take the velocity of the center-of-mass of the two fluid 
locally as v , that is, v Q = U, we have 

AC/ b = U> + - Ul = i( 7 -C7 + -7+C/-) = J± ^(V + - V_), 

7 7 



and 



\AU h \ < \U' + - U'_ 
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We can write 



U± = — ( 7 ±pC/±7TPAC/ b ), (130) 

P7 

K = ^V + ^^AU b . (131) 

en n 2 



In general, since A/i < 1, the condition of the neglect of the current momentum is 

eTl |,r t t I 1^1 



Pe<— , |E/+-E/_|<^. (132) 



This is because 



J 


= l-^l = 




ne 




en 



^AU» 

n z 



< P±\ V \ + ^±^|AC/ b | < ^|t/| + \U + - U-\ « \U\. 

ne n l ne 



The influence of the unique nature [3] is estimated by comparison between A/i j (ne) 
and in the second term, the right-hand side of equation (|4"8l) . As the above consideration, 
if 

Pe «g, |tf + -U_|«M, (133) 

the unique nature[3]is negligible. In the pair plasma case (A/i = 0), the Hall effect disappears 
without any restriction. 

The distinct property [21 the time-derivative term of the Ohm's law, is evaluated by the 

ratio 

^j^—UJ : WF* V ~ J : B. (134) 

cr [ney ct {ne) z 

The ratio should be evaluated in a more realistic case because it depends on the characteristic 
time scale r. An example on black hole magnetospheres of AGNs is discussed in the next 
section. With respect to the unique property H] of the thermal electromotive force, the 
influence is estimated by the ratio 

where I is the characteristic scale length. As for the unique nature the term is negligible in 
the electron-ion plasma case in a short time scale, because of the energy exchange between 
the electron fluid and ion fluid. However, we cannot ignore it in the pair plasma case. In 
this paper, we do not discuss the details of the last two distinct properties. 
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7. Validity of generalized and standard RMHD equations 

In this section, we show the applicability of the generalized RMHD equations and evalu- 
ate the significance of the pair and electron-ion plasma distinct properties in a global, astro- 
physical situation. The standard RMHD equations are applicable only when the conditions 
of the generalized RMHD are satisfied and the unique properties of pair and electron-ion 
plasmas can be ignored (see section [6]). 



7.1. Break-down condition of generalized RMHD equations 



We consider the break-down condition of the generalized RMHD equations fH6T) -(H£ 
with the Maxwell equations Q, (|1J) on the basis of the two-fluid approximation and the 
average procedure fl5}-((HJ). To derive the generalized RMHD equations from the one- fluid 
equations, we used the assumptions M = U U U U = — 1 and AhJ <C h. It is noted that the 
one-fluid equations (12"TT) - (I23"I) and the two-fluid equations ([I])-© are consistent within the 
averaging procedure (|oTl-([5|). The condition M = — 1 is necessary so that has meaning 
of the four-velocity, and that means the averaged Lorentz factor of the fluid of positively 
charged particles and the other fluid of negatively charged particles is non-relativistic in 
the center-of-mass frame because M = U U U U = UjU u ' = Uq'U ' = — ( / y 1 ) 2 = — 1. Here, 
the prime indicates a variable observed by the local center-of-mass frame of the two fluids. 
This condition is consistent with the assumption of the current-density independency of the 
resistivity (see Appendix A). When the relative velocity of the two fluids is relativistic, the 
frictional force of the two fluid is not proportional to the relative velocity, and then the 
resistivity depends on the current density. This condition restricts the net charge density 
and current density as follows. Using equations © and dHJ, we get 



P 1 



(n+ - n) 2 + -J U J V 



by some algebraic calculations. When M 



U V U V w 1, we find 



n_ 



< i. 



(136) 



Using this condition and J V J U = —pi + \ J\ 2 , we find 



n_ 



then we have 



e z 



e 2 m_i_m_ 



\J\ 2 « pl 



[en\ 



P 



(137) 
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This condition shows the current must be much smaller than en/ \J~jl or p e . 

For an electron- ion plasma, if <C m + /m_ and n + pa n_, we have the condition 7' « 1. 
This is because 

Tfl 

l' ~ 7+ H -7- ~ 1, 

m + 

since 

|E/' + | = ^|t/'i<— 7' «1. 

m + n + m + 

For a pair plasma, we need 7^ ~ 1 to keep the condition of 7' pa 1 because 7^ ~ 7^ and 
7' ~ (7+ + 7-)/2 ~ 7± if n + « n_. 

With respect to the condition Ah' <C ft,', in the case of the non-relativistic pressure 
limit p± <C p±, we find 

( \ 1 \ m 1 m_ 
2pm < — + . 

\n + n_ J n + n_ 

This condition requires n + pa n_. It is noted that this does not mean the charge neutrality 
because net charge density is given by p e = e(n + 7 + — n_7_) 7^ e(n + — nJ). We call the 
condition n + pa n_ the "proper charge neutrality". In the relativistic pressure case, we find 
this condition yields 

2 f h+.m_ h^m+\ h + h_ 

< -4 + 



m \ ni n_ I n± n_ 



and then using specific enthalpy (I65I) . which is approximated by equation (1661) . we have 

J.* C ^ „ ±.h ( P -) . (138) 



n + \P+J n- VP-, 
This condition is satisfied when n + ~ n_ and p+/p+ ~ p_/ p_ except for a special case. 

In the case of an electron-ion plasma with proper charge neutrality, the condition 
p+/ p+ pa p-/ p- is implausible when the ion temperature T + is not much larger than the elec- 
tron temperature T_ because it requires T + = p + /n + pa (m+/m-)m-P-/p- = (m+/ra_)T_ 
and m + ^> m_. Then, the application condition of the generalized RMHD equations is 
p± <C p±. In a case of a pair plasma with proper charge neutrality, the application con- 
dition is p+/p+ ~ p-l ' P- or p ± <C p ± . In general (both cases of the electron- ion and pair 
plasmas), the generalized RMHD equations are applicable to the plasmas with the proper 
charge neutrality n + pa n_, the non-relativistic thermal energy p± <C p±, and the non- 
relativistic (non-superrelativistic for electron-ion plasmas) relative velocity of the two fluids 
U' + < U'_ <C m + /m_ if m + > m_. 

The premise conditions of the generalized RMHD equations are summarized as follows. 
The generalized RMHD approximation breaks down if all of the following conditions are not 
satisfied. 
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1. The local relative velocity of the two fluids is non-relativistic, 7' ~ 1. This condition 
comes from the condition U U U U = — 1. This condition restricts the current density as 
shown by equation (I137p . In an electron- ion plasma, it is satisfied when 7^ <C m + /m_ 
and n + w n_, while in a pair plasma it requires 7^ ~ 1. 

2. The condition of the proper charge density n + « n + is required to hold the condition 
A/;,"'" <C . It is noted that this condition does not mean the charge neutrality p e ~ 0. 
In the non-relativistic pressure case (p± < the condition Ah^ <C b) is satisfied 
only with this condition. 

3. In the relativistic pressure case, we require the condition shown by equation (11381) . In 
the electron-ion plasma case, this condition yields p± <C p± except for the special case 
such as T + « (m + /m_)T_. In the pair plasma case, it requires T + « T_ or p± <C p±. 



7.2. Unique nature significance of electron-ion/pair plasmas in AGN jet 

forming region 

We investigate whether unique properties of the electron-ion and pair plasmas shown 
by their generalized RMHD equations are significant in a jet forming region around an AGN 
black hole, which has not been done by previous works with standard RMHD equations. 
Here, we assume that in such a region global magnetic field crosses an accretion disk around 
a black hole and the plasma around the blac k hole has zero resistivity, which is a typical as- 



sumption for jet formation models (see FigJH) ( Blandfo rd fc Payne 19821; IShibata fc Uchida 1986 ; 



Kudoh et al. 19981 : iKoide et al. 19981 : iKoide et al. 19991 : iKoide et al. 20001 ) . The disk twists 



the magnetic field lines rapidly near the black hole to increase the magnetic pressure. The 
magnetic pressure (and tension) blows off the plasma near the disk. The outflow of the 
plasma is collimated by the magnetic tension to form a jet. We here make a rough esti- 
mation of the distinct properties of the electron-ion and pair plasmas through the process 
around the black hole. For simplicity, we ignore the general and special relativistic effects 
such as the frame-dragging effect of the rotating black hole and the inertia effect of the 
magnetic field. The former effects are discussed in section [SJ In this section, we use the 
MKSA units. 

As shown in Fig. [1], in such a situation, current is induced along the rotation axis 
of the disk, and the return current is formed. The generalized RMHD equations suggest 
that the current carries the momentum to influence the jet formation. We evaluate the 
influence as follows. We assume the velocity of the disk rotation is the Keplerian velocity 
vk = a/ GMbr/t, where G is the gravitational constant, Mbh is the mass of the black 
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hole, and r is the distance from the center of the black hole to the jet forming region. The 
azimuthal component of the twisted magnetic field B& is estimated by 



B l 1 2 

" -Vk, 



2p p 2 

where p is the mass density of the coronal plasma around the black hole. On the other 
hand, the Ampere's law yields poJ ~ B^/r. As shown in section [6] (equation (11281) ). we 
evaluate the influence of current momentum and the Hall effect by the ratio [/z|«7|/(ne)]/|L/| 
and pp e / \nerf). The former ratio is estimated as 



1 J c e m cp 



— —p . — -p 

V ne r y ne z ru p 

In the pair plasma case, the Hall effect disappears. In the electron-ion plasma case, it is 
negligible when 

c 



Ap pa < 1. (139) 



The critical particle number density, with which the current momentum or the Hall effect is 
significant, n cr it, near the black hole is estimated by 



—J 6 ° m max(/i, Ap,) ~ 1, 
rs V ^crite^ 



and using max(p, Ap) < 1, we have 



"""it ~ 3^2 



r s e 



° 2 C ° m {max(/i, Ap)} 2 < f / W = 3.2 x 10 6 ( [m^] . 

1 VP, P) S _ (2GMbh) 2 e 2 \ AfgH / L J 



Here, we set r = rg where rs = 2GMbh/c 2 is the Schwarzschild radius. In the case of an 
AGN, we calculate the critical number density as n cr i t ~ 3 x 10~ 2 [m -3 ] when Mbh ~ 10 8 M Q . 
This critical number density n cr j t is much less than the actual value of the plasma around 
the black hole and this evaluation shows insignificance of the current momentum and Hall 
effect. 



We evaluate the inertia effect of the current in the Ohm's law (jUj). As shown in section 
|6] (equation (11341) ). the influence is estimated by the ratio 

p mn J p c ( c \ t> K 



cr n 2 e 2 B Tu p ru p \ruj p 



where r is the characteristic time scale of the jet formation r ~ r/v^- This shows if the 
current momentum and the Hall effect are negligible, we can neglect the current inertia effect 
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because of equation (11391) and the condition vk < c. As the results, the distinct properties 
of the pair and electron-ion plasmas in the jet forming region near the black hole are not 
significant globally in the AGN case. Inversely, the condition that the current momentum is 
significant is 



The Hall effect is not negligible when 



r < — fi. 



c 

r < — A// 



in the electron-ion plasma case, while it disappears in the pair plasma case. The inertia of 
the current is significant in the case that 



1 fCOJp 



Here, we neglect the thermal electromotive force and the thermalized energy exchange be- 
tween the two fluids. These conditions can be satisfied only in the small-scale phenomena. 
The global plasmas of black hole magnetospheres may be influ enced by these disti nct effects 
through the magnetic reconnection around the black holes (IKoide fc Arai 20081 ). This is 
discussed in section M briefly. 



8. Generalized GRMHD equations for electron-ion/pair plasmas around black 

holes 



Around a foot-point of a jet near a black hole shown in Fig. 1, the general relativistic 
effect becomes significant because of strong gravity of the black hole. In this section, we 
discuss the general relativistic effects of the plasma dynamics near black holes, while in the 
previous sections we investigated the RMHD equations in the special relativistic framework. 
To treat the plasmas near black holes appropriately, we have to use the general relativistic 
equations. We get the equations by transforming the generalized RMHD equations (1461) - 
with the Maxwell equations (j3J), (jlj) to the covariant form of the general relativity 



V v (pU v ) = 0, 
V 



h ( UW + ^J"J U 



2fiAfih 



(140) 
(141) 



Q 



^V(A/ip - Ap) + (U u - ) F''„ //[./" ! Q( 1 OH '"i 

2q y q 



(142) 
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V„ = 0, 



(143) 
(144) 



where V M is the covariant derivative (IMisner et al. 19701 ; IWeinberg 1972T ). Using the equa- 
tion derived by the Maxwell equations 

1 



{V v F IUT )F m = -y^v yr r nX ) 
and equation (11441) . the equation of motion (11411) becomes 



where 



JlflV 



P9 



fj,i/ 



+ h (uw 



V„T^ = 0, 



+ J»J V ) + F\F ua - -y» v {F 



(145) 



kX f kX ) 



This equation corresponds to the equation of motion in the standard GRMHD, for example, 
equation (A2) in Appendix A of Koide et al. (2006). The newly additional term is only that 
of the current momentum density pJ^J u /q 2 . The Ohm's law becomes 



V, 



2fiAfih 



Q 



-V^(Afip -Ap)+[U 



A/i 



J v F% - //[J" + Q(l + Q)W 



2q \ q J 

when we ignore the non-uniformity and time variation of q. As shown in section [3j H = 
{rjq) 2 j (ph) should be smaller than several (almost unity) because of causality. Equation (11461) 
with H < 1 is called the equation of the "generalized general-relativistic Ohm's law" and 
equations QUOD , ffMTD . ffT46l) . ffT43D . fTM with H < 1 are called the "generalized GRMHD 
equations" . Equation (I146p is drastically different from the ideal MHD condition U U F' 1 U = 0, 
which is used customarily in present GRMHD simulations. To show the plasma dynamics 
governed by the generalized GRMHD equations, we discuss them using the 3+1 formalism 
shown by Koide, Kudoh, & Shibata (2006) (see Appendix A of the paper; definition of 
variables therein). The variable observed by the "zero angular momentum observer (ZAMO)" 
frame is denoted by a hat. The energy density observed by the ZAMO frame is 



e + 7p = T 



hj 2 • ^ P \-v • 



B 2 E 2 



where e+jp is the total energy density. This shows "electric charge" has the rest mass energy, 
which is disregarded in the standard GRMHD. The three components of the momentum 
density, 



P l = T w = h\ -iU l + ^pj 

qz 



[E x B) 
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show that the current has the momentum. The stress-tensor, 

f-y = P 5v + h (u'w + ^J'J^j + + 613 ~ ^ j ~ ^ v 

indicates the current also carries the momentum. Summarizing the above results, we find 
that the difference of equation of motion between the generalized GRMHD and standard 
GRMHD equations is only with respect to the inertia and momentum of the current, which 
are found in the special relativistic MHD equations. 

We next examine the Ohm's law. Defining 
we have another formalism of the general relativistic generalized Ohm's law 

v,T o 7 m = ^(a^ - Ap) + (ir - *£jr) f\ - w + q(i + e)cn (147) 

The i-th component of equation (I147P yields the 3+1 form 

dt \q 2 / hh 2 h 3 ^ ° 3 [ hj 1 ohm + P q 2 



fih 1 da 



q 2 hi dx i 



pi 



+ - £ + E { G *r pt - G *&) + Tai^ W ~ AP) (M8) 



+ 1 [E i + (V x B)i] - ^ [ Pe E l + ( J x B)J - v [.r + Q(l + B)U% 

where h 2 , h 2 ,, h\ are spacial diagonal elements of the metrics, /3 2 , /3 3 are the shift vectors, 
a is the lapse function, is the shear of the frame-dragging = (l/hj)[d(a{3 1 ) /dx^], 
fk m = Z^GijfiL - Cj.V'D- and G t] = Here, p\ = (q 2 7 >)T o h ° m = 

2p e ( 7 - App c /(en)) and J ]i = (q 2 /ph)f^ m = jj { + p c U l - (2Ap/q)pJ i are on the order of 
the charge density and the current density, respectively, and a is regarded as gravitational 
potential. With respect to the term —(ph/q 2 )(l/hi)(da/dx l )pl, it clearly shows the gravity 
influences the Ohm's law and behaves like electromotive force as Khanna (1998) pointed out. 
Near the black hole where the gravitational acceleration is very strong, such "gravitational 
electromotive force" becomes significant when the charge separation of the plasma occurs. 
This effect is also found in the Newtonian two-fluid model when we consider the gravitation 
term in equation of motion of the two fluids. The forth and fifth terms of the right hand 
side of equation f!148j) 

7ti , \^„.mP h ( 



q* * — ' q z 

3 3 
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is also regarded as the effective electromotive force. In the electromotive force, the shear of 
the frame dragging around black holes, (fih/q 2 )aji and (fjih/q^ceftGik, corresponds to the 
electric resistivity. This effective resistivity due to the frame-dragging shear is a pure general 
relati vistic effect, which may cause the dynamo effect and magnetic reconnection nea r black 
holes (Khanna fc Camenzind 1994; Khanna fc Camenzind 1996: Koide fc Arai 2008). 



9. Summary and discussion 

In this paper, we derived the one-fluid equations of the two-component plasma from 
the relativistic two-fluid equations. On the basis of the new one-fluid equations and the 
dispersion relation analysis of the electromagnetic wave, we proposed a set of generalized 
RMHD equations ([IB" ]) — (j4~8"l) with Maxwell equations ([3]) and (jl]) which are applicable to 
both pair and electron-ion plasmas. The generalized RMHD equations yielded the RMHD 
equations for the electron-ion plasma (equations (14"9"]) - (15T]) ) and the pair plasma (equations 
( 155|) - (l57j) ). The comparison between these sets of equations and the standard equations 
( I52p - (|54p clearly show the distinct properties of the electron-ion and pair plasmas. We 
investigated the linear modes propagating in the pair and electron-ion plasmas using the 
generalized RMHD equations. In the analysis, we found the effect of polarity rotation of 
Alfven wave in the electron-ion plasma, while no rotation in the pair plasma. This effect is 
similar to the Faraday rotation of the electromagnetic wave propagating along magnetic field 
lines in a plasma, but not identical. We also found other plausible unique properties of the 
electron-ion and pair plasmas, such as the plasma oscillation, fast wave, and electromagnetic 
wave, which also show the validity of the generalized RMHD equations. With respect to the 
non-linear effects, we evaluated the distinct properties of the pair and electron-ion plasmas 
with the significance of each distinct term of the generalized RMHD equations and with the 
special case of black hole magnetospheres where jets are formed. It confirmed the validity 
of the standard RMHD equations fl52|) - fl54^) . ([3]), (j4j) on the global phenomena in black hole 
magnetospheres of AGNs. We also discussed the generalized GRMHD equations f)140p - 
( 11441) and revealed the unique effects of the plasmas around black holes. As indicated by 
Khanna (1998), we confirmed the appearance of the "gravitational electromotive force" due 
to the strong gravity of a black hole when the plasma has charge separation. We newly 
found the effective resistivity due to shear of frame dragging around a rotating black hole. 
To investigate the non-linear phenomena qualitatively, numerical method is required and 
should be developed in the near future. 

According to the results in subsection 17.11 if the relative velocity of the two fluids 
and the thermal energy of the fluids are non-relativistic and their proper particle number 
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densities are approximately equal, 7' ~ 1, p± <C p±, n + pa n_, the generalized RMHD 
equations are applicable. Otherwise, the premise condition (A/J <C ftJ) of the generalized 
RMHD equations require p+/p+ ~ P-/P-, that is, T + = (m + /m_)T_. This condition is 
implausible for the electron-ion plasma when the ion temperature is not much larger than 
the electron temperature. In the pair plasma case, it yields p + /n + = T + pa T_ = 
The special condition of the electron-ion plasma, T + = (m + /m_)T_ may be possible for 
optically thin Advection Dominated Accreti on Flows (ADAFs) in the electr on-ion plasma of 



the accretion disk around black holes (e.g.. iKato Fukue fc Mineshige 19981 . see Table 10.1). 



If these premise conditions are not satisfied, the generalized RMHD equations, which are 
reducible to the standard RMHD equations, are not applicable to the plasma. That is, a 
relativistically hot plasma (p± > p±) can not be treated by RMHD if p + pa p^{m + /mJ) and 
n + pa n_ are not satisfied. In such a case, the generalized RMHD approximation breaks down. 
Then, the generalized RMHD equations should be replaced by more appropriate equations 
for the plasma with relativistic internal energy (relativistic high temperature, relativistically 
fast relative velocity of the two fluids). The standard RMHD equations, which have been 
utilized indiscriminately in astrophysics, are applicable to the astrophysical objects when the 
premise of the generalized RMHD equations are satisfied and the unique properties of the 
electron- ion/pair plasma are negligible. These conditions indicated in this paper depend on 
the averaging procedure of the variables ffTOl) . The averaging of the variables fflUj) enables 
us to derive the one-fluid equations (T2TT) — (T23T) from the two-fluid equations ([I])-© relatively 
easily. However, we have to check the validity of the averaging method. To confirm it, the 
Vlasov-Boltzmann equation of the plasma should be utilized. This task is beyond the scope 
of this paper, while it is interesting and important. It should be investigated in the near 
future. 

As shown in subsection 17.21 distinct properties of the electron-ion and pair plasmas 
do not influence the global dynamics of the plasmas in black hole magnetospheres of AGNs 
because they are large-scale RMHD phenomena. On the other hand, local phenomena such as 
magnetic reconnection may change the influence to the global dynamics drastically, because 
they are possible to change the magnetic configuration globally, although the reconnection 
region is much smaller than the global scale of magnetospheres. In such a small region, 
the distinct properties of the electron- ion and pair plasmas may become significant. For 
example, the inertia effect of the current plays a role to cause the magnetic diffusion instead 
of the electric resistivity. This distinct magnetic reconnection may be important in the pair 
plasma, which is thought to locate at the coronal regions near black holes of AGNs and 
GRBs. 



With respect to the corona between a black hole and its disk, in generally speaking, the 
plasmas would be composed of electrons, ions, and positrons, while the amount of positrons in 
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the plasmas around black holes has not been confirmed observationally yet. This is regarded 
as a mixture of electron-ion plasma and pair plasma. The simplest treatment of the mixed 
plasma is given by m_ = m e and m + = m c + (m ; — m e )£, where m- x is the mass of the ion and 
£ = rii/n e - is the relative proportion of ions (n ; and n e - are the particle number densities 
of ions and electrons, respectively). Here, we assumed the local (proper) charge neutrality, 
n e - = rii + n c +, where n e + is the number density of the positrons. In this treatment, the 
relative proportion of ions £ should be given by other condition. For example, we may be 
able to trace the plasma element where £ is constant. However, when the relative velocity 
between the positron fluid and the ion fluid is significant, £ on a certain plasma element 
changes. In such a case, the treatment becomes difficult. The advanced formalism of the 
mixed plasma is a forthcoming and important subject in the relativistic plasma physics. 
Additionally, other effects such as pair creation/annihilation, radiation, atomic processes, 
and more general EoS should be considered. Furthermore, as we mentioned above, for the 
relativistically hot plasma and relativistically strong current plasma, the premise of the 
generalized RMHD equations breaks down. In plasmas at black hole magnetospheres, it is 
not surprising that the temperature and current density are relativistically high and strong. 
To deal with such relativistic plasmas, we have to develop a new scheme beyond the present 
generalized RMHD equations. 

I am grateful to Mika Koide, David L. Meier, Masaaki Takahashi, and Takahiro Kudoh 
for their helpful comments on this paper. This work was supported in part by the Sci- 
ence Research Fund of the Japanese Ministry of Education, Culture, Sports, Science, and 
Technology. 



A. Frictional four-force density of two fluids 

We assume a relative velocity of two fluids is non- relativistic so that M = U„U U — — 1. 
Then, when we observe the fluids from the local center-of-mass frame S', the fluids can be 
treated as non-relativistic fluids. The four-force density of the friction from the negatively 
charged fluid to the positively charged fluids is 

/;' = -m + a + _v T n + n_y + y_(Vi! - V!!), (Al) 

where <t + _ is the cross-section of the collision between the positive particle and negative 
particle, v T is the average relative velocity between the two fluids including the thermal 
velocity. Similarly, the frictional four-force density of the positively charged fluid to the 
negatively charged fluid is 

ft' = -m_a_ + v T n + n_i + i_{V* - V»). (A2) 
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According to action-reaction principle, we have 

m + o" + _ = m_o"_ + . (A3) 
Then we define the following variables, 

A = r/e 2 = m + cr + _v r = m_a- + v r . (A4) 
In the center-of-mass frame S', we have 

Ul = 0. (A5) 
When we consider the Lorentz transformation x M = b^^x"', we find 

b% = —. (A6) 

7 

We calculate the frictional four-force density by 



7 7 



j» _ Q_ uf . 

M 



where M = U V U V = - 7 /2 and Q = U U J U = ~i J Q ' . We also note that 

/|V = - 2 ^' /2 (Q 2 " MW) < 0, 



where W = J U J U and 



/iw = - 2 Apm r, 2 (Q 2 -Mw). 

m z n_y_ye z 



We consider the thermal energy gain rate of the positive charged fluid, 

f+ = i+n+ + i _ n _ n +^+ ~ f+ v +* t 1 - 9 +) - ("A v +* ) 



7+n + + 7^n_ 
Ap 9 + m^[pM - (m + /e)Q} 2 - 9^m + [pM + (m_/e)Q] 



-(Q 2 -MW), 



me 2 (2pM + ^±Q) (pM - =*Q) (pM + =^Q) 

where 6 I + , are the redistribution coefficient of the thermalized energy to the positively 
charged fluid and negative one with the equipartition principle, respectively. When 9 = 9 + = 
9-, we obtain 

If = < 0, (A7) 

me 
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where 



= 



9 (Q 2 - MW) [(m+ - m_) {-(pM) 2 + 



4pMQm + m- 

e 



eQ 



(2p M + ^±g) (pM - =*Q) (pM + Z±Q) 



Using 0, we finally obtain the frictional four-force density, 



Ap 



ft = -V^M 



me 



J"- ^(i + e)c/" 



Because i? M = /+, we obtain the Ohm's law, 
1 



2ne 
1 



2ne 



(2ne) 5 



^(Ap - Ap) + [U v -J v F" v - r]V^M 

2ne \ ne 



Q 



J" - — Ci + eW" 



(A8) 



(A9) 



(A10) 
• (AH) 



Here, we recognize rj is the resistivity. 
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Fig. 1. — Schematic picture of a relativistic jet forming region in a black hole magnetosphere 
of an AGN. The disk is thought to be mainly composed of a electron-ion plasma, while the 
corona near the black hole and the core of the jet would be mainly of pair plasma. The 
coronae near the disk and outer shell of the jet probably consist of mixed plasma with ions, 
electrons, and positrons. Magnetic field lines are twisted by the accretion disk rotating 
around the black hole, and the jet is driven by the magnetic pressure. In this case, the 
current and return current are formed around the jet, which is indicated by the arrows of 
the current density J. 



